Measurement of Azimuthal Asymmetries in Inclusive Production 
of Hadron Pairs in e^e Annihilation at Belle 
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The Collins effect connects transverse quark spin with a measurable azimuthal dependence in the 
yield of hadronic fragments around the quark's momentum vector. Using two different reconstruction 
methods we find evidence of statistically significant azimuthal asymmetries for charged pion pairs in 
e'^e~ annihilation at a center-of-mass energy of 10.52 GeV, which can be attributed to a transverse 
polarization of the primordial quarks. The measurement was performed using a sample of 79 million 
hadronic events collected with the Belle detector. 

PACS numbers: 13.88.+e,13.66.-a,14.65.-q,14.20.-c 



The relationship between quark spin and the properties 
of hadrons is still poorly understood. In the fragmenta- 
tion of a quark into hadrons, Collins P| has proposed 
a relation between transverse quark spin and the final 
state azimuthal distribution of hadrons around the origi- 
nal quark momentum direction. Azimuthal asymmetries 
of hadron yields have recently been reported in deep in- 
elastic scattering off transversely polarized hydrogen 
and deuteron |2j targets: these are beheved to be due to 
the product of Collins fragmentation functions and trans- 
verse quark spin distribution functions in the nucleon. 
However, these so-called transversity distributions, which 
contribute to the nucleon transverse spin, have not been 
independently measured. Similarly, no measurement of 
Collins fragmentation functions has yet been reported. 

In this Letter we present a measurement of azimuthal 
asymmetries in hadron-hadron correlations for inclusive 
charged dihadron production e~^e~ — > hhX, which we 
interpret as a direct measure of the Collins effect. Ex- 
ploratory studies for this measurement were reported 
in Q. This analysis was performed using a sample of 79 
million events (29fb~^) collected at a center-of-mass sys- 
tem (CMS) energy 60 MeV below the T(4S') resonance 
with the Belle detector at the KEKB asymmetric- 
energy e+e" storage rings [^. For systematic checks, 
Monte Carlo (MC) simulated events generated by the 
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QQ and JETSET packages and processed with a full 
GEANT-based simulation of the Belle detector were 
used. 

The Collins effect occurs in the fragmentation of a 
quark q with transverse spin Sq and 3-momentum k into 
an unpolarized hadron h with transverse momentum Ph_L 
with respect to the original quark direction. The corre- 
sponding number density is defined as 

D,^r (z, P,^) = Dfiz, Pl^)+Hi\z, Pl^) ^ • , 

zMh 

(1) 

where z = is the fractional energy of the hadron rela- 
tive to half of the CMS energy Q. The first term describes 
the spin averaged fragmentation function (FF), and the 
second, containing the Collins function H^'^{z, P^^), de- 
pends on the spin of the quark and produces an asymme- 
try as it changes sign if Sq flips. The vector product leads 
to a sin((/)) modulation, where (j) is the azimuthal angle 
between the plane spanned by the hadron and quark mo- 
menta, and the plane spanned by the quark momentum 
and the incoming leptons. Experimentally the quark di- 
rection is approximated by the thrust axis n. 

In hadron production in e+e^ qq events, the Collins 
effect can be observed when the fragments of the quark 
and anti-quark are considered simultaneously. Combin- 
ing two hadrons from different hemispheres in jet-like 
events, with azimuthal angles 0i and 02 as defined in 
Fig. n (note that all angular vairablels as well as n are de- 
fined in the CMS), produces a cos{<j>i + (j>2) modulation of 
the di-hadron yield. A MC comparison of thrust axis cal- 
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Thrust axis n 



FIG. 1: Definition of tiie azimutiial angles of tiie two hadrons. 
In each case, is tiie angle between the plane spanned by the 
lepton momenta and the thrust axis n, and the plane spanned 
by n and the hadron transverse momentum Phi±- 



culations using reconstructed and generated tracks shows 
an average angular deviation between the two of 75 mrad, 
with a spread with root mean square of 74 mrad. This 
smearing of the reconstructed axis leads to a reduction in 
the measured azimuthal asymmetry, as discussed below. 

Two experimental methods are used to measure az- 
imuthal asymmetries. The first method {M12) gives rise 
to the cos(0i + 02) modulation in the di-hadron yields. 
The yield is recorded as a function of the hadron angle 
sum (l)i+(l>2, N\2 = Ni2{4)i-\-(j)2), and divided by the aver- 
age yield to obtain the normahzed rate R12 := Ni2{cj)i + 
4'2)/{Ni2), parametrized by i?i2 = ai2 cos((^i -I- (^2) + 612- 
Here, 012 is a function of the first moment (iJ-j^'^''^') of 
the Collins function lldl 



ai2{0,Zi,Z2) 



(^l)^l 



^2 



1 



COS'' 



DKzi)Dl{z2) 



(2) 



where 9 is the angle between the incoming lepton axis 
and the thrust axis. An alternative method (Mq) does 
not rely on knowledge of the thrust axis: yields are mea- 
sured as a function of (jjQ, the angle between the plane 
spanned by the momentum vector of the first hadron and 
the lepton momenta, and the plane defined by the two 
hadron momenta. The corresponding normalized rate 
Ro = Nn(2(t>n) / (No) is a function of cos(2<^o), and (fol- 
lowing can be parametrized as aq cos(20o) + with 



00(612,2:1,22) 



(2l)i//(z2)/MlAf2) 



1 



cos" 



D\{z,)D\{z2) 



(3) 

/ denotes convolution over the transverse hadron mo- 
menta. Ml and M2 are the masses of the two hadrons, zi 
and Z2 are their fractional energies and 62 is the angle be- 
tween the beam axis and the second hadron momentum. 
The sin 02 dependence reflects the probability of finding 
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FIG. 2: Top: Unlike(U)-sign and like(L)-sign pion pair nor- 
malized rate _Ro vs. 2<^o in the bin 2:1(2:2) € [0.5, 0.7], 2:2(21) £ 
[0.3,0.5]. Bottom: Pion pair double ratio Rq /Rq vs. 2(^o in 
the same bin. The solid and slashed lines show the results of 
the fit described in the text. 



the two initial quarks with transverse spin. Di(z) and 
Hf"^ denote fragmentation functions for anti-quarks. 

To reduce hard gluon radiation, a two-jet-like topology 
is enforced by requiring a thrust value T > 0.8, calcu- 
lated from all charged and neutral particles with momen- 
tum exceeding 0.1 GeV/c. The following selection criteria 
were imposed on the charged pions used in the analysis 
methods M12 and Mq: (1) Tracks are required to origi- 
nate from the collision vertex, and to lie in a fiducial re- 
gion —0.6 < cos(6'iab) < 0.9, where 6'iab is the polar angle 



in the laboratory frame. (2) A likelihood ratio is used to 
separate pions from kaons P : C{n)/[C{K)+C{'K)] > 0.7. 
MC studies show that less than 10% of pairs have at least 
one particle misidentified. (3) We require zi, Z2 > 0.2, to 
reduce decay contributions to the pion yields. In ad- 
dition we require the visible energy in the detector to 
exceed 7 GeV. (4a) The tracks must lie in opposite jet- 
hemispheres: (Phi ■n)(P/i2 -n) < 0. (4b) Qt is the trans- 
verse momentum of the virtual photon from the e~^e~ an- 
nihilation in the rest frame of the hadron pair jlll ]. We 
require Qt < 3.5 GeV/c, which removes contributions 
from hadrons assigned to the wrong hemisphere. 

The analysis is performed in (zi, Z2) bins with bound- 
aries at Zi = 0.2, 0.3, 0.5, 0.7 and 1.0, where complemen- 
tary off-diagonal bins (21,22) and (22,21) are combined. 
In each (21, 22) bin, normalized rates i?i2 and Rq are eval- 
uated in 8 bins of constant width in the angles 4>i + (f>2 
and 200 respectively, and fitted with the functional form 
introduced above. Results in the lowest (21,22) bin are 
shown in Fig. [3 In both methods the constant term (612 
or bo) is found to be consistent with unity for all bins. 

In addition to their sensitivity to the Collins effect, i?i2 
and Rq have contributions from instrumental effects and 
QCD radiative processes: these are charge independent. 
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FIG. 3: Values of Ao and A\2 as functions of z, corrected for 
the contribution of charm events. The lower scales show the 
boundaries of the bins in zi and Z2; see the text. The shaded 
band shows the size of the systematic errors. 



and cancel in the double ratio of normalized rates for 
unlike-sign(U) over like-sign(L) pion pairs, 

Bl _ N^{I3^)/{N^) 

a = 0,12, /3o = 20o, /3i2 = 01+02- (4) 

In linear approximation the double ratio is proportional 
to a combination of the favored and disfavored fragmen- 
tation functions: omitting the transverse momentum de- 
pendence, 



R^/R^ = l + cos(20o 



sm 



1 + cos'^9 



(5) 



an analogous expression can be given for RY2/ Ri2- Fa- 
vored fragmentation processes [e.g. D(°'^) are transitions 
in which a valence quark of the hadron is of the same fla- 
vor as the initial quark, for example u,d tt"*"; the corre- 
sponding unfavored process is u,d —^ 7r~ . Following our 
analysis of the normalized rates we parametrize the dou- 
ble ratios as R^/R^ ~ cos(/3a) -I- Ba with a = 0, 12, 
f3o = 200 and P12 = 0i -I- 02 ; the parameters Aa , Ba are 
determined from fits in each (zi, Z2)-bin. 

An example is shown in Fig. |21 The parameters Aq 
and A12, related to the Collins fragmentation functions 



(see Eq. are shown in Fig. |21 and listed in Table 
Significant non-zero values are observed, especially at 
high z. The observed increase with z is predicted in sev- 
eral models |0, Q| • The weighted averages over all 
(2:1, Z2) -bins, corrected for the charm contribution (see 
below), are found to be Aq = (3.06 ± 0.65 ± 0.55)% and 
A12 = (4.26 ± 0.78 ± O.68)%0. 

By measuring asymmetries using double ratios, accep- 
tance effects cancel, while the contribution of gluon radia- 
tion cancels only to first order. The size of any additional 
contribution was estimated by subtracting like-sign from 
unlike-sign pair rates, in which case gluon effects should 
cancel exactly and only experimental effects should re- 
main. The differences, 0.04% on average for method Mq 
and 0.03% for Afi2, are assigned as systematic errors. 
Double ratios were formed using MC events, which do 
not include the Collins effect but take into account gluon 
radiation and detector effects: the parameters Aq and 
A12 were found to be consistent with zero. The statis- 
tical errors of these fits, 0.45% and 0.56% respectively, 
are included in the systematic uncertainty. Double ra- 
tios of positively-charged over negatively-charged pairs 
were found to be consistent with unity, and limit the 
possibility of charge dependent detector effects; again, 
the precisions of the fits (0.26% and 0.21%) are included 
as systematic errors. An additional test was performed 
by taking pion pairs from jet hemispheres in different 
events: no asymmetry was found. 

In addition to fitting the double ratios with a cos0 
modulation, higher harmonics (sin 20, cos 40) were intro- 
duced. The change in the results, 0.03% on average for 
Aq and 0.01% for A12, was included in the systematic un- 
certainty. In method M12, the single hadron yield around 
the thrust axis was also studied. Although the Collins ef- 
fect leads to a sinusoidal modulation, this should average 
to zero in the absence of a specified quark spin. The 
sin 01^2 modulation was found to be consistent with zero. 

As a consistency check, we also measure double ratios 
from an event sample with a reversed thrust selection, 
T < 0.8: asymmetries in this sample are reduced for 
Qt < 3.5 GeV/c (see Fig.0J. For the low-thrust sample 
there is no clear two-jet topology, and, thus, the Collins 
effect should be suppressed while any radiative or ac- 
ceptance effects remain present. A fit to a constant for 
Qt < 3.5 GeV/c finds (0.4 ±0.3)% for both Aq and A12. 
The asymmetries for Qt > 3.5 GeV/c and T < 0.8 were 
found to be due to events with low visible energy as well 
as background from interactions in the detector material. 

The uncertainty due to particle misidentification was 
estimated by applying tighter selection criteria. The dif- 
ference in the double ratio with respect to the default 
selection is included in the systematic error, 0.12% for 
method Mq and 0.28% for M12. 

To test the reconstruction of azimuthal asymmetries, 
Monte Carlo samples were reweighted in cos(20o) and 
cos(0i -1-02), producing generated moments of 5% and 
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TABLE L Aq and A12 values obtained from fits to pion double ratios as a function of z. The errors shown are statistical and 
systematic. 





Zi ^ Z2 




Ao 


Ar2 


[0.2,0.3] 




[0.2,0.3] 


(1.68 ±3.10 ±0.54)% 


(1.08 ±3.53 ±0.67)% 


[0.2,0.3] 




[0.3,0.5] 


(0.77 ± 1.46 ± 0.54)% 


(2.74 ± 1.76 ± 0.67)% 


[0.2,0.3] 
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(3.36 ± 1.49 ± 0.54)% 


(3.60 ± 1.80 ± 0.67)% 
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(4.10 ±1.95 ±0.54)% 


(2.84 ±2.37 ±0.67)% 
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(2.71 ± 1.82 ± 0.54)% 


(1.34 ±2.21 ±0.67)% 


[0.3,0.5] 
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(5.19 ±1.56 ±0.54)% 
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[0.3,0.5] 




[0.7, 1.0] 
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(6.37 ±2.37 ±0.67)% 
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(5.24 ± 4.26 ± 0.62)% 
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(21.84 ± 8.34 ± 1.59)% 
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FIG. 4: Values of Aq and A12 as a function of Qt (not 
corrected for charm background) for the data samples with 
T > 0.8 (triangles) and T < 0.8 (squares). 



10% for unlike-sign pairs, and 0% and —5% for like-sign 
pairs. The reconstructed azimuthal asymmetries were 
consistent with generated values for method Mq, but 
lower by 17.6 ± 1.1% on average for method M12, which 
depends on the thrust axis of the reconstructed event. 
(Note that the thrust axis represents the direction of the 
outgoing quarks only on average.) This dilution was cor- 
rected by rescaling A12 with a factor 1.210 ± 0.014. 

Both methods depend on the assumption that the 
electron and positron beams are unpolarized. This was 
tested by studying the angular distribution of e+e^ 
/Lt+/z~ events: no significant azimuthal asymmetries are 
observed and thus no systematic error assigned to it. 
We correct for the contribution from charm decays us- 
ing measured asymmetries in events where a D* meson 
has been reconstructed, and the e+e~ cc event frac- 
tion determined from MC (23%). The corresponding un- 
certainty is included in the statistical errors in Table 
The fraction of selected events due to the e+e" t^t" 
process is small (1.7%), and the asymmetries obtained in 
a tau enhanced data sample have low statistical signifi- 
cance: Ao = (0.17 ± 0.30)% and A12 = (0.74 ± 0.30)%. 



This contribution is added to the systematic error. All 
systematic uncertainties were added in quadrature. Cor- 
relations among individual angular and (zi, 2:2)-bins were 
tested using a large number of MC samples, by compar- 
ing uncertainties returned by the fits with the expected 
values. We find the statistical error to be underestimated 
by 14%; the final uncertainty is increased by the corre- 
sponding factor. 

In summary we have performed a measurement of the 
azimuthal asymmetry in the inclusive production of pion 
pairs as a function of the fractional energy z of pion pairs. 
In the double ratio of asymmetries from unlike-sign and 
like-sign pairs, possible contributions of gluon radiation 
and detector effects cancel and the observed asymmetry 
can be attributed to the Collins effect. This asymmetry 
in e^e~ annihilation is the first direct evidence for this 
effect. 
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